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Abstract. We prove C 2 ' a regularity of sufficiently flat free boundaries, for 
the thin one-phase problem in which the free boundary occurs on a lower 
dimensional subspace. This problem appears also as a model of a one-phase 
free boundary problem in the context of the fractional Laplacian (—A) 1 / 2 . 



1. Introduction 

Let g(x, s) be a continuous non-negative function in the ball B\ C R ,l+1 = 
R™ x R, which vanishes on a subset of R n x {0} and it is even in the s variable. We 
consider the following free boundary problem 

f Ag = 0, in B+(g) := B 1 \ {(x,0) : g(x,Q) = 0}, 

(1-1) { dq 

(-^j=l, on F(g) := ^{x e Bi : g(x, 0) > 0} n B u 

where 

(1.2) —{x ):=hm + — , x G F{g) 

with v(xq) the normal to F(g) at xo pointing toward {x : g(x, 0) > 0} and B r C R™ 
the n-dimensional ball of radius r (centered at 0) . 

If F(g) is C 2 then it can be shown (see Section 7) that any function g which is 
harmonic in B^(g) has an asymptotic expansion at a point xq € F{g), 

g(x, s) = a(x )U((x - x ) ■ v(xq), s) + o(\x - x | 1/2 + s 1/2 ). 

Here U{t, s) is the real part of y/z which in the polar coordinates 

t = rcosd, s = rsin6>, r > 0, — tt < 8 <tt, 

is given by 

n 

(1.3) U(t,s) = r 1/2 cos-. 

Then, the limit in (ll.2[) represents the coefficient a(xo) in the expansion above 
(which justifies our notation) 

J^(zo) = a{x ) 

and our free boundary condition requires that a = 1 on F{g). 

Solutions to our free boundary problem (jl.l[) are critical points to the energy 
functional 

E(g)-= f \\7g\ 2 dxds + ^n n ({g>0}nB 1 ). 
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If the second term is replaced by H n+1 ({g > 0}), we obtain the classical one-phase 
free boundary problem (see for example |ACj .) In our case the free boundary occurs 
on the lower dimensional subspace M. n x {0} and for this reason we refer to (jl.ip as 
to the thin one-phase free boundary problem. 

This free boundary problem was first considered by Caffarelli, Roquejoffre and 
Sire [CafRS] as a model of a one-phase Bernoulli type free boundary problem in the 
context of the fractional Laplacian. It is relevant in applications when turbulence 
or long-range interactions are present, for example in flame propagation and also 
in the propagation of surfaces of discontinuities. For further information on this 
model see |CafRS| and the references therein. 

In this paper we are interested in the question of regularity for the free boundary 
F(g). Concerning this issue the authors of }CafRS] proved that in dimension n = 2, 
a Lipschitz free boundary is C 1 . In [DR], the first author and Roquejoffre showed 
that in any dimension if the free boundary F(g) is sufficiently flat then it is C l ' a . 

This paper is the first of a series of papers, which investigate the regularity 
of F(g) and in particular the question of whether Lipschitz free boundaries are 
smooth. This basic question was answered positively in the case of minimal surfaces 
by De Giorgi fDG] and by Caffarelli jClj for the standard one-phase free boundary 
problem. 

Our strategy to obtain the regularity of Lipschitz free boundaries is to use a 
Weiss-type monotonicity formula jW] combined with flatness results and ad hoc 
Schauder type estimates near the free boundary. To implement this method we 
need to obtain first C 2a estimates for flat free boundaries, which we achieve in this 
paper. Unlike the case of minimal surfaces and of the standard one-phase problem, 
C 2 ' a estimates do not seem to follow easily from C 1,a . It appears that C 2 ' a is the 
critical regularity needed to obtain C°° smoothness of the free boundary, as well as 
the regularity needed to implement our blow-up analysis. 

The following is the main result of this paper (see Section 2 for the precise 
definition of viscosity solution to (jl.ljl ). 

Theorem 1.1. There exists e > small depending only on n, such that if g is a 
viscosity solution to Ijl.ljl satisfying 

{x E B x : x n < -e} C {x E B x : g(x, 0) = 0} C {x e B x : x n < e}, 

then F(g) is a C 2 ' a graph in Bi for every a G (0, 1) with C 2,a norm bounded by a 
constant depending on a and n. 

The proof of Theorem 11.11 follows the lines of the flatness theorem in [DR , 
which is inspired by the regularity theory developed by the second author in [Sj. 
In this case the proof is more technical since we need to approximate the free 
boundary quadratically. To do so, we introduce a family of approximate solutions 
Vs,a,b which have the same role as quadratic polynomials in the regularity theory 
of elliptic equations. Such family will be used also in a subsequent paper to obtain 
boundary Schauder type estimates for solutions to our problem. 

In the last section of this paper we also prove some useful general facts about 
viscosity solutions g to our free boundary problem (|1.1[) . such as C 1 / 2 -optimal 
regularity, asymptotic expansion near regular points of the free boundary and com- 
pactness. 

The paper is organized as follows. In Section 2 we recall notation, definitions 
and some basic results from [DRj . including the linearized problem associated to 
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(jl.lj) . Section 3 is devoted to the construction of the quadratic approximate so- 
lutions Vs, a ,b- In Section 4 we prove a Harnack type inequality for solutions to 
In Section 5 we establish the improvement of flatness result via a compact- 
ness argument which makes crucial use of the Harnack inequality of Section 4. Our 
argument reduces the problem to studying the regularity of solutions to the lin- 
earized problem. This is pursued in Section 6. We conclude the paper with Section 
7 where we provide some general facts about viscosity solutions to (jl.ljl . 

2. Definitions and basic lemmas 

In this section we recall notation, definitions and some necessary results from 
(DR]. 

2.1. Basic facts. Throughout the paper, constants which depend only on the di- 
mension n will be called universal. In general, small constants will be denoted by 
c, Ci and large constants by C, Ci and they may change from line to line in the body 
of the proofs. 

A point X e K™ +1 will be denoted by X — (x, s) € R™ x R, and sometimes 
x = (x' , x n ) with x' = {x\, . . . , x n -x). 

A ball in R n+1 with radius r and center X is denoted by B r (X) and for simplicity 
B r = B r (0). Also B r denotes the n-dimensional ball B r n {s = 0}. 

Let v <E C(Bi) be a non-negative function. We associate to v the following sets: 

B+(v) := B 1 \ {(x,0) : v(x,0) = 0} C 

B+(v) :=B 1 + (t))nBiCM n ; 

F(v) :=d*nB£(v)nBi CM™. 

Often subsets of R n are embedded in R™ +1 , as it will be clear from the context. 
We consider the thin one-phase free boundary problem 



A<7 = 0, in B+(g), 
(2-1) { do 



where 



— x := hm T , X = {x , G F(g). 

aU t^o+ V< 



Here v(xq) denotes the unit normal to F{g) 1 the free boundary of <?, at xq pointing 
toward B^{g). 

We now recall the notion of viscosity solutions to (|2.1[) . introduced in [PR] . 

Definition 2.1. Given g,v continuous, we say that v touches g by below (resp. 
above) at X S B x if g(X ) = v(X ), and 

g(X) > v(X) (resp. g(X) < v(X)) in a neighborhood O of X . 

If this inequality is strict in O \ {Xq}, we say that v touches g strictly by below 
(resp. above). 

Definition 2.2. We say that v g C(B\) is a (strict) comparison subsolution to 
(|2.ip if v is a non-negative function in B\ which is even with respect to s = and 
it satisfies 
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(i) v is C 2 and Av > in Bf(v); 

(ii) F(v) is C 2 and if x G F(v) we have 

u(a;o + ti/(xo), 0) = o^ieo)^ + o(v / i), as < ^ + , 

with 

a(x Q ) > 1, 

where v{x$) denotes the unit normal at xq to F (v) pointing toward (v ); 

(hi) Either v is not harmonic in B±(v) or a(a;o) > 1 at all xq £ F(v). 

Similarly one can dehne a (strict) comparison supersolution. 

Definition 2.3. We say that g is a viscosity solution to (|2.f jl if g is a continuous 
non-negative function in _Bi which is even with respect to s = and it satisfies 

(i) Ag = in B+(g); 

(ii) Any (strict) comparison subsolution (resp. supersolution) cannot touch g 
by below (resp. by above) at a point X — (xq, 0) € F(g). 



Remark 2.4. We remark that if g is a viscosity solution to (|2.1j) in B\, then 

SaPO = A- 1 / 2 5 (AX) ! X e By 
is a viscosity solution to (|2.ip in Si. 



Finally, we state for completeness the boundary Harnack inequality which will be 
often used throughout the paper. This version follows from the boundary Harnack 
inequality proved in [CFMS . 

Theorem 2.5 (Boundary Harnack Inequality). Let v be harmonic in B^~ (v) and let 

F(v) be a Lipschitz graph in the e n -direction (pointing towards the positive phase) 
with £ F(v). If w is harmonic in B^"(w) = B^~(v), then 

— < C— -e„ in B 3/i , 
v v 2 ' 

with C depending only on n and on the Lipschitz constant of F(v). 
2.2. The function g. Here and henceforth we denote by P the half-hyperplanc 
P :={X £ R n+1 :x n <0 lS = 0} 

and by 

L := {X £ R n+1 : x n = 0, s = 0}. 

Also, throughout the paper we call U(X) := U{x n ,s), where U is the function 
defined in flO). _ 

Let j be a continuous non-negative function in B p . As in [PR] , we define the 
multivalued map g which associate to each X £ M ,l+1 \ P the set g(X) C M via the 
formula 

(2.2) U(X)=g{X-we n ), Vw £ g(X). 

We write g(X) to denote any of the values in this set. 

This change of variables has the same role as the partial Hodograph transform for 
the standard one-phase problem. Our free boundary problem becomes a problem 
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with fixed boundary for g, and the limiting values of g on L give the free boundary 
of g as a graph in the e n direction. 

Recall that if g satisfies the e-flatness assumption 

(2.3) U(X - ee n ) < g(X) < U(X + ee n ) in B p , for e > 

then g(X) ^ for X e B p - e \ P and |s(^0| < e, hence we can associate to j a 
possibly multi-valued function <? defined at least on -B p _ c \ P and taking values in 
[— e, e] which satisfies 

(2.4) U(X) = g(X - g(X)e n ). 

Moreover if g is strictly monotone in the e n -direction in Bj{g), then g is single- 
valued. 

We recall the following lemmas from [PR] . 

Lemma 2.6. Let g,v be non-negative continuous functions in B\ with v strictly 
increasing in the e n -direction in Bt(v). Assume that g and v satisfy the flatness 
assumption (|2.3[) in B\ for e > small. If 

v < g in B\, 

then 

v < g on B X -e \ P. 

Viceversa, if 

v < g on B a \ P, 

for some < a < X — e, then 

v < g on B a - C . 

Lemma 2.7. Let g,v be respectively a solution and a subsolution to (|2.1|) in Bi, 
with v strictly increasing in the e n -direction in B^iv). Assume that g and v satisfy 
the flatness assumption (|2.3p in B^ for e > small. If, 

(2.5) v + o-<g in(B 3/2 \B 1/2 )\P, 
for some a > 0, then 

(2.6) v + a < g in B 3/2 \ P. 

Finally, given a Lipschitz function (j) defined on B\(X), with values in [—1,1], 
then for all e > small there exists a unique function ip e defined at least on B\^ e (X) 
such that 

(2.7) U{X) = ^{X-e4>{X)e n ), X e B X (X), 
that is 

ip € = 6<ft. 

Moreover such function (p e is increasing in the e„-direction. 

If g satisfies the flatness assumption (|2.3[) in B\ and <fi 1S as above then (say 
A < 1/4, X e Bi /2 ,) 

(2.8) &<9 mB x (X)\P^ip e <g in B X ^{X). 

The following Proposition will be used in the compactness argument for the proof 
of the improvement of flatness in Section 6. 
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Proposition 2.8. Let <j> be a smooth function in B\(X) C K™ +1 \ P. Define (for 
e > small) the function (p e as above by 

(2.9) U(X) = <pt{X-vf>{X)e n ). 
Then, 

(2.10) A^ e = eA{U n <j>) + 0(e 2 ), in B X/2 (X) 
with the function in 0(e 2 ) depending on \\<fi\\c 5 an d A. 

Proof. For notional simplicity we drop the subindex e in the definition of ip e . From 
formula (|2.9p and Taylor's theorem, we have that 

(2.11) U{X) = ^{X)-c Vn {X)4>{X)+cH{X), mB x/2 (X) 
with ||^'||c 3 (b a/2 (a > )) < C and C depending on ||0||c 5 an( i ^- Thus, 

U n (X) = <p n (X) + 0(e). 
Combining this formula for ip n (X) and (|2.11|) we obtain 

U(X) = y{X) - eU n (X)cf>(X) + 0(e 2 ). 
Hence, using that U is harmonic, 

= AU(X) = A<p(X) - A{eU n cl>)(X) + 0(e 2 ), 
as desired. □ 

We remark that in fact the function in 0(e 2 ) only depends on A if we choose e 
small enough depending on ||<^||c5. 

2.3. The linearized problem. We recall here the linearized problem associated 
to (|2.ip . Here and later U n denotes the a;„-derivative of the function U . Recall that 

P:={X e R n+1 :x n <0,s = 0}, L :— {X e : x n = 0, s = 0}. 

Given h G C{B{] and X = (x' , 0, 0) £ Bi fl L, we call 

IV7 u\tv\ T ^(^O) s ) — ^(^Oi 0, 0) 2 2,2 

V r /i (A ) := hm , r = x + s . 

(x„,s)->(0,0) r 

Once the change of unknowns (|2.2j) has been done, the linearized problem associated 
to (HU) is 



(2.12) 



A(U n h) = 0, in B l \ P, 
\\7 r h\ = 0, on Bi n L. 



Definition 2.9. We say that h is a solution to (|2.12j) if h <E C(£?i), ft- is even with 
respect to {s = 0} and it satisfies 

(i) A{U n h) = in Si \ P; 

(ii) h cannot be touched by below (resp. by above) at any X — (x' Q ,0,0) £ 
Bi n L, by a continuous function cf> which satisfy 

4>{X) = cb(X ) + a(X ) ■ (x' - x' Q ) + b(X )r + 0(\x' - x' \ 2 + r 3 / 2 ), 

with b(X ) > (resp. b(x ) < 0). 

In Section 6, we will prove a quadratic expansion for solutions to the linearized 
problem which yields the following corollary. 
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Corollary 2.10. Let h be a solution to (|2.12|) such that \h\ < 1. Given any 
a £ (0, 1), there exists 770 depending on a , such that h satisfies 

\h{X) - (MO) + Co • x' + \{x') T M Q x' - |V - b rx n )\ < m B Vo , 

with r 2 = x\ + s 2 ; for some a ,b Q £ R, £ £ R n_1 ,Mo G S(n-l)x{n-l) with 

\£ \, \ao\, \bo\, \\M \\ <C, C universal 

and 

ao + bo — trMo = 0. 



3. A FAMILY OF FUNCTIONS. 

In this section we introduce a family of functions Vs, a ,b which approximate our 
solution quadratically. These functions will be often used as comparison subsolu- 
tions/supersolutions. We establish here some of their basic properties, including 
their behavior under the change of coordinates V — > V (see Proposition 13. 5jl . 

We start by presenting some basic properties of the solution U defined in the 
introduction. Recall that 

U(t,s) ^p^cos^, 

where 

t = pcos(3, s — p sin f3, p > 0, — tt < (3 < it. 
We will use the following properties of the function U : 

(i) AU = 0, U > in M" +1 \ P. 

(ii) u t = \p~ 1/2 cos f = —U and U t > in M™ +1 \ P. 

lp 

Since Ut is positive harmonic in IR 2 \{(i, 0), t < 0}, homogenous of degree —1/2 
and vanishes continuously on {(£,0), t < 0} one can see from boundary Harnack 
inequality (or by direct computation) that values of Ut at nearby points with the 
same second coordinate are comparable in diadic rings. Precisely we have 

(3.1) 777r4^ C if \ti-t2\<ht 2 ,s)\. 

Ut(t2,s) 2 

Next we introduce the family V5 j0j f,. For any 0,6 £ I we define the following 
family of (two-dimensional) functions (given in polar coordinates (p, /?)) 

(3-2) v a , b (t, s) := (1 + ^p + h -t)p 112 cos |, 



that is 



!>„, 6 (t, s) = (1 + ^p + h -t)U{t, s) = U(t, s) + oip 1 / 2 ). 



Given a surface S — {x n = h(x')} C M™, we call Ts,x the 2D plane passing 
through X — (x, s) and perpendicular to S, that is the plane containing X and 
generated by the s-direction and the normal direction from (x, 0) to S. 

We define the family of functions 

(3.3) V s , a ,b{X):=v a>b {t,s), X = (x,s), 
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with t = pcos/3, s = psin/3 respectively the first and second coordinate of X in the 
plane Vs.x- In other words, t is the signed distance from x to S (positive above iS 
in the ir n -direction.) 
If 

S:={x n =^(x') T Mx' + e-x'}, 
for some M e x e we use tne notation 

(3.4) V M ,e>, a ,b(X) := V s<a< b{X). 
This will be the case throughout most of the paper. 

Definition 3.1. For S > small, we define the following classes of functions 
V S := {V M ,?, a , b : \\M\\,\?\,\a\,\b\<6}, 

and 

V? - {VM,$>, a ,b e V s : a + b- trM = 0}. 

Most of the times we will work with functions in the class Vs , since we deal with 
the flat case. Notice that if we rescale V = Vmx' ,«•*> that is 

V X {X) = \- 1 ' 2 V{\X), X G Bi, 

then it easily follows from our definition that 

In the next proposition we provide a condition for a function V £ Vs to be a 
subsolution/supersolution. 

Proposition 3.2. Let V = VM,£,',a,b € V$, with 5 < So universal. There exists a 
universal constant Cq > such that if 

(3.5) a + b- trM > C 5 2 
then V is a comparison subsolution to (12 . 1|) in B2 ■ 

Proof. Clearly from our formula for v a & the function V satisfies the free boundary 
condition of Definition 12.21 with a(xo) = 1. We need to check that AV(X) > at 
allXe B+(V). 

Since that V(X) depends only on (t,s) and 

A x t = —k(x) 

where k{x) is the sum of the principal curvatures of the parallel surface to S (in 
W 1 ) passing through x, we compute that 

(3.6) AV(X) = A(t,s)v a ,b - (d t v a j } ) K {x). 
From our formula for v a ,b, using polar coordinates we get that 

(3.7) A (t)S) u a , b = i(a + typ- 1 ' 2 cos | = (a + b)U t . 
Also, since p < 2, 

(3.8) |Sle««,6 ~ C7*| < (|o| + IH)p 1/2 cos | < 8<5C/ f . 
Finally we use that «i(a;) the principal curvatures at x are given by, 
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□ 



where x* is the projection of x onto S. Since |£'|, ||M|| < 6 we obtain that 

Ma*)| < OS, \k(x*) - trM\ < CS 3 
for C universal, which in view of (|3.9[) give 

(3.10) \k(x) - trM\ < CS 2 . 
From (|3J]l combined with (|37f|) . (pTgjl and (|3~T0l) we get that 

(3.11) |AV(X)- (a + 6-trM)C7 t | < icV 2 [/ 4 

for a Co universal. It follows that if 

a + b - trM > C 5 2 

then AV(X) > as desired. 

Next, we estimate V n and AV outside a small cone with axis L. 
Proposition 3.3. Let V = VM,£',a,b € Vs with S < 5q universal, then 

(3.12) c <Y^<c, inB 2 \(PU{\(x n ,s)\<106\x'\}). 
IfV e V 5 ° then 

(3.13) \AV{X)\ < CS 2 U n {X) inB 2 \(PU{\(x n ,s)\ < 106\x'\}). 
Proof. From our formula 

dt 

V n {X) = d t v a , b {t,s) — 

OXji 

where t represents the signed distance from x to S. Since V ' x t is the unit vector at 
x that has the direction of the normal from x to S, it makes an angle of order <5 
with respect to e n . Hence since 

dt _ t 

OX n 

we get 

(3.14) 1 > > 1 - CS 2 

ox n 

and we obtain 

d t V a .b(t,s) > V n (X) > -d t V atb (t,s). 

From p.8j) we see that dtv a .b ~ Ut and we obtain that 

(3.15) 2d t U(t, s) > V n {X) > ^d t U(t, s). 

Thus to obtain our claim we need to replace t with x n in the inequality above. 

Since in B 2 \ x \ the surface iS is in a 4S\x\ neighborhood of x n = we find that 
I* — %n\ < 4<J|x|. If X belongs to the domain in (|3.12[) then 

Kan,*) I > 8<5|a;| > 2\t - x n \ 

and we obtain from (13. ip 
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which together with (|3.15[) gives the desired conclusion (|3.12p . 

Now (|3.13[) follows immediately. Indeed by formula p.ll|) we have that 

\AV(X)\ <CS 2 U t (t,s), 

which combined with (13.161) gives the desired bound. □ 

Remark 3.4. We remark that if V G V$, then the rescaling V\(X) = X~ 1 ^ 2 V(XX) 
with A < 1, satisfies 

(3.17) c < ^ < C, \AV X (X)\ < C6 2 U n (X), 

in the dilation ball of factor 1/ A 

B 2/x \(PU{\(x n ,s)\ < 10S\x'\}). 

Indeed 

AV X (X) = X 3 ^ 2 AV(XX), U n (X) = X^ 2 U n (XX), (V x ) n (X) = X^ 2 V n (XX). 

Now we study the behavior of V G Vg under the transformation V — > V. This 
will be quite useful in the rest of the paper. 

Proposition 3.5. Let V — VM,£',a,b G Vs, with 6 < 5q universal. Then V is 
strictly monotone increasing in the e n -direction in £?^"(V). Moreover, V satisfies 
the following estimate in B 2 \ P 

\V(X) - lv {X)\ < C^ 2 , MX) = \r 2 + brx n - \{x') T Mx' - ? ■ x' 

with r = \fx% l + s 2 and C\ a universal constant. 
Proof. First we show that v a .b satisfies 

(3.18) U(t + 7a , b - CS 2 , s) < v a>b (t, s) < U(t + 7a , fc + CS 2 , s) 

where p 2 = t 2 + s 2 and j a .b is the following expression depending on t and s: 

. . a 2 
7a,fe(*,sJ := -p +bpt. 

Indeed since (see properties of U listed at the beginning of this action) 

\U tt \ <Cp- x \J t 

we have that if \p\ < p/2 then 

\U(t + p, s) - (U{t, s) + pU t (t, s))\ < p 2 \U tt (t', s)\ < Cp 2 p' 1 U t (t, s), 
where in the last inequality we used (I3.1[) . Thus, since Ut = U/(2p), 

(1 + f + c£)U(t, s) > U(t + M) s) > (1 + f - C^)U(t, s), 
2p p A 2p p z 

Choosing 

p = p ± 4C^ 
P 

we obtain that 

~ 2 ~ ~ 2 

U(t + p + 4C^-,s) > (i + JL)U(t,a) > U(t + p-AC^,s), 
p 2p p 
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provided that \p/p\ < c, with c sufficiently small. Since 

a b . 

we can apply the inequality above with 

a i 

V=2 P + btp > 

hence \p\/p < CS and obtain the claim. 

When t is the signed distance from x to the surface S we have 

t = on S := {x n = h(x') := -x' T Mx' + £' ■ x'} 

and by (|3~T4l) 

1>^>1-CT 2 in Bx, 
da; n 

thus, by integrating this inequality on the segment (x', h(x')), [x 1 ,x n ) we get 

\t-{x n -h{x'))\ <CS 2 . 

Since in B\, the surface S and x n — are within distance 5 from each other we 
have \t — x n \ < CS and hence 

|7o,b(*) s) - y a ,b(%n, s)\ < \\Vv a ,b\\L°° \t - x n \ < CS 2 . 

From the last two inequalities we have that 

\(t + la,b{t, a)) - {x n + 1V (X))\< CS 2 , 

with 

lv {X) = y a , b (x n , s) - ^x' T Mx' - g ■ x'. 

Using this fact and ()3.18[) (and the monotonicity of U in the e n direction) we 
obtain 

U(X + ( 1V (X) - CS 2 )e n ) < V(X) < U(X + ( lv (X) + CS 2 )e n ), 

and the estimate for V is proved. 

Finally, we remark that the monotonicity of V follows from (|3.15[) . □ 

Remark 3.6. Notice that from the last inequality in the proof above, we obtain that 
if V € Vg, then V satisfies the 45-flatness assumption in B\ (see also (|2.3t ): 

U(X - 4Se n ) < V(X) < U(X + ASe n ). 

This could be also checked easily directly from the definition of V. 

We conclude this section with by comparing the functions V corresponding to 
two nearby surfaces. 

Lemma 3.7. Let Si,i = 1,2 be surfaces with curvature bounded by 2. Let 
Vi = V Suai M> N,N<2, i = 1,2. 

Assume that, 

Si n B 2(7 = {x n = hi(x')}, a < c 
with hi Lipschitz graphs, hi(0) = 0, |V/lj| < 1 and c universal. Lf 
\ai - a 2 \, \h -b 2 \ < e, \\hi - h 2 \\L^ < tcr 2 , 
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for some small e < c, then 

Vi(X) < V 2 (X + Cea 2 e n ) in B a . 

Proof. After a rescaling of factor 1/cr, we need to prove our lemma for a = 1 and 
with the curvature of Si, dj, bi and e smaller than c universal. 
First we prove that for < A < 1, 

«oi,ii(t,a) <v a2 , b2 (t + Ce\ 2 ,s), \<p = \(t,a)\ < 2A. 

By (|3.15|) , dtVa.b is proportional to <9tC7 in the disk of radius 2. Since on the 
segment with endpoints (t, s) and (t + CeX 2 , s) all the values of d t U are comparable 
(see (|3.ip ) we obtain (using 2pC/ t = U) 

v a2 M( t + CeX 2 , s) > v a2 M 2 {t, s) + Ce\ 2 U t {t, s) 
>f/(l + fp+|* + Ce^) 

> U 0l ,6i (*,*)> 

and our claim is proved. 

Since v a2 i, 2 is increasing in the first coordinate, we obtain that 

Uai.biCM) < Wa 2 ,b 2 (i + Ce, s), |(t,s)| < 1. 

On the other hand, from the hypotheses on hi we see that in B\ 

h + Ce< i 2 , 

where <2 is the distance to 52 — C'ee n , for some C large depending on the C above. 
Hence in B\ we have 

Vi(X) = w Ql!bl (ti,s) < u 02 ,b 2 (*i + Ce,s) < w a2 ,6 2 (i 2 ,s) = V 2 (X + C'ee n ). 

□ 

4. Harnack Inequality 

In this section we state and prove a Harnack type inequality for solutions to our 
free boundary problem (|2.1[) . This will allow us to obtain some compactness of 
flat solutions after the transformation g —> g (see Corollary I4.2j) which is a crucial 
ingredient in Theorem ll.il 

Theorem 4.1 (Harnack inequality). There exist e > small and C > large 
universal, such that if g solves (|2.1I) and it satisfies 

(4.1) V{X + a e n ) < g{X) < V{X + b e n ) in B P (X*) C B u 
with V — VM,£',a,b & Vg, and 

CS 2 < b -^-^ < l, 
P 

with \oq\, |6o| < 1) then 

(4.2) V{X + a 1 e n )<g{X)<V(X + b 1 e n ) in B^ p {X*), 
with 

a < ax < bi < b , 6 X - a x = (1 - fj)(b - a ), 
for a small universal constant fj G (0, 1/2). 
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In the particular case when V = U, this statement was proved in }DR| . Our proof 
follows the same lines as the one in [PR] but it requires a more careful analysis since 
the function V is no longer a precise solution. 

From this statement we get the desired corollary to be used in the proof of our 
main result. Precisely, assume g satisfies (|4.1I) in B\ with do = — e, bo = e for some 
small e <C e, and S such that CS 2 < e. Notice that from Remark l3.6[ the functions 
V and g are (4(5 + e)-flat in B\. 

Then at any point X* £ -B1/2 we can apply Harnack inequality repeatedly for a 
sequence of radii p m — ^fj m and obtain 

V(X + a rn e n ) < g(X) < V(X + b m e n ) in B^ m (X*), 

with 

(4.3) b m - a m = (b - oo)(l - = 2e(l - ^) m , 

for all m's, m > 1 such that 

(4-4) 4 ^ < e. 

This implies that for all such m's, the function g satisfies 

(4.5) V + a m <g< V + b m , in Si fim _ 45 _ £ (X*) \ P, 

with a m ,b m as in (|4.3j) . Define the following (possibly multivalued) function 

(4.6) g 6 , v pO := 9(X)-V(X) ^ ^ g Bl _ 44 _ £ \ Pj 
and notice that 

In view of (j4~5|) we then get that in Bi^ (X*) \ P 

(4.7) osc g e , v < 2(1 - f]) m , 
provided that 

(4.8) 4S + e<e 1/2 <fj m /4. 

If e < e fj 2m " for some nonnegative integer mo then our inequalities above (|4.4I) . 
(|4.8p and hence also (|4.5[) hold for all m < mo. We thus obtain the following 
corollary. 

Corollary 4.2. Let g solve (|2.ip and satisfy for e < e 

V{X - ee n ) < g(X) <V{X + ee n ) in B u 

with 

V = V M , M € Vl CS 2 <e, 
for e, C > universal constants. If 

e <efj 2mo , 

for some nonnegative integer m (with fj > small universal), then the function 
g t y defined in (|4.6p satisfies 

a e (X)<g €tV (X)<b e (X), mB 1/2 \P 

with 

b t -a t <2(\~n) m °, 
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and a € ,b e having a modulus of continuity bounded by the Holder function at^ for 
a, /3 depending only on fj. 

The proof of Harnack inequality will follow from the Proposition below. 

Proposition 4.3. There exist e, S > and C > universal, such that if g solves 
(|2.1[) and it satisfies 

(4.9) V(X-ee n ) < g(X) < V(X + ee n ) in B u for < e < e 
with 

V = V M ,s', a ,b eV,°, 5 < 5, C5 2 < e, 

then either 

g(X)<V(X + (l-r ) )ee n ) in B v , 

or 

g > V(X - (1 - r?)ee„) mB„ 
for a small universal constant r\ € (0, 1). 

First we show that if g > V and they separate of order e at one point, then they 
separate also of order e away from a neighborhood of L — {x n = 0, s = 0}. This 
follows from the boundary Harnack inequality. Below are the details. 

Lemma 4.4. If g solves (I2.1[) and it satisfies 

(4.10) g(X) > V(X - ee n ) in B u 

(4.11) g(X)>V(X) at some X E Bi(-e n ), 

with V = VM,t', a ,b € Vj;, C<5 2 < e for C > universal, then 

(4.12) <?(X) > y(X - (1 - r)ee„) C, 
with 

C := {(x',x n ,s) : - < \(x n ,s)\ < -, \x'\ < -}, d: 



2 - lx " - 2" 1 - 2 J ' 8^^! 
and t a small universal constant t £ (0, 1). 

Proof. We have 

V(X - (1 - r)ee„) = V^(X - ee„) + reK(Jf + Aee„), 

for some A with |A| < 1. Hence by (|3. 1[) . ()3. 12|) for e small enough and X E C 

V(X - (1 - r)ee„) < F(X - ee„) + CTeU n (X + Aee„) 

< V(X-ee n ) + CtTeU n (X). 

Thus, if h(X) := g(X) — V(X — ee n ) we need to show that 

(4.13) h > aeU n , in C, 

and then choose t = c\/C\. 

To obtain ()4.13|) . notice that by a similar computation as the one above in view 
of (|4~TT|) and (|3~Tj) . (f3~i"2"l) we get that for e small enough 

(4.14) h(X) > V(X) - V(X - ee n ) > cU n (X)e > c 2 e. 
Also, by (|4.10p we have 

h>0 in Si. 
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Finally, by (pTl3)l 

\Ah\ < CS 2 U n < C 2 S 2 , in C\P, 

where C DD C is the d/4-neighborhood of C. 

Thus in view of (|4.14p and Harnack inequality we have that (for C large enough) 

(4.15) h > c 2 e - CS 2 > c 3 e, in B 1/8 (~e n ). 
Denote by 

D:=C\(B 1/8 (±e n )UP) 

and let q\ , q 2 satisfy in D 

(4.16) Agi = 0, Aq 2 = -1 
with boundary conditions respectively 

qi = on dC U P, 31 = 1 on d-Bi/g(-e„) 

and 

(72 = on 9D. 

By boundary Harnack inequality, q\ is comparable to the distance function s in a 
neighborhood of P(~\C CC C. Since q 2 is Lipschitz continuous in a neighborhood of 
P n C, we then obtain 

(4.17) q x > c 4 q 2 in C \ B 1/8 (~e n ), 

with C4 > universal. By the maximum principle, 

h > q := c 3 eqi - C 2 S 2 q 2 in D, 

since h > q on and Ah < Aq in D. Hence, by (|4.17|) we get that (for C large 
enough) 

h > eygi > c 5 eU n , in C \ B 1/s (-e n ), 

where in the last inequality we used that (by boundary Harnack inequality) q\ and 
U n are comparable. This inequality together with (I4.15|) gives the desired claim 

flfn]). 

□ 



We are now ready to present the proof of Proposition 14.3 
Proof of Proposition \4-3\ Assume that 



(4.18) g(X)-V(X)>0, X= l -e n . 

Then in view of assumption (|4.9j) from Lemma 14.41 after the change of variables 
g — >• g we get that 

(4.19) g{X) > V(X) +re-e in C \ P 
with 

C :={{x',x n ,s):d<\{x n ,s)\<-. \x'\<-}, d ' 



4" '"2 J ' 87^1 
Denote by 

W{X) := V M +^eI,Z',a,b+2ce(X) G V 5+e , 
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with c small to be made precise later. Then in view of Proposition 13.51 we have 

(4.20) - 20,(6 + e) 2 <(V-W) + ce(2rx n - 1 \x'\ 2 ) < 2d(S + e) 2 . 

2(n — 1) 

First we choose c small depending on t such that 

V> W- T -e, 

where we used that CS 2 < e < 1 with C > C(r) and e small enough. Then, if C is 
sufficiently large depending on c, 

(4.21) V>W + r*e, on{\(x n ,s)\<d,\x'\ = ^}\P, 

for some r* > small, say t* < r/2. These combined with (14. 19|) give 

(4.22) g(X)>W(X)+T*e-e in (C U {\(x n , s)| < d, \x'\ = i}) \ P. 

Moreover, if (7 is large enough we get that W satisfies (|3.5p and hence 14 7 is a 
subsolution. Thus from Lemma 12.71 and the inequality above we conclude that 



(4.23) g(X)> W(X) + r*e-e in {\(x n , a)\ < d, \x'\ < ^} \ P. 

Finally, from (I4.20[) we see that there is a small neighborhood around the origin 
B v C {K^nj s )| < d, \x'\ < |} ijj small universal depending on the constants above, 
r] < r*/2) such that 

W > V - T — e, in B 2v \ P. 
Hence, from (|4.23|) we conclude that 

g > V + m - e in B 2v \ P, 

for some small universal constant r), and the lemma is proved after the change of 
variable g — > g. □ 

We conclude this section with the proof of Theorem 14.11 



Proof of Theorem After a translation of the origin we may assume that we 
satisfy our flatness hypothesis (|4.ip in B p (X*) c B% with 

(x*)' = 0, a + & = 0, VeV° s . 

We dilate the picture by a factor of 2/p and work with the rescalings 

9 P (X) = (f r 1/2 g(^x), V P (X) = fy-V'V^X), 

which are defined in a ball of radius 2 included in B^/ p . Notice that, if V G 
then V p £ V° s . 

After dropping the subindex p for simplicity of notation, we may assume that 
the flatness condition (gTj holds in some ball B 2 (X*) C W l+1 , with V G V° 5 , 

a = -c, b Q = e, (x*)' = 

and 

C7) 2 < e < e. 

We need to prove the conclusion (|4.2I) in a ball B2fj(X*). 
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We distinguish three cases depending on whether X* is close to L, close to P, 
or far from P. 

In Case 2 and Case 3 we will use the following properties from Remark 13.41 
(4.24)c<^<C, \AV\<C6 2 U n in B 2 (X*) \ (P U {\(x n , s)\ < 206\x'\}). 



Below r\ is the universal constant from Proposition 14.31 
Case 1. \X*\ < rj/A. 

In this case, since P>i C B 2 (X*) we follow under the assumptions of Proposition 
Hence we can conclude that for any fj < rj/A in B 2 fj(X*) C B v either 

g(X) < V(X + (1 - r,)ee n ), 

or 

g(X) > V(X - (1 - r,)ee n ), 
and our conclusion is satisfied for all fj < r] / A. 

Case 2. \X*\ > r]/A, and B^_{X*) n P = 0. 

In this case, if e is small enough then it follows from (|4.24p that the function 
h{X) :=g(X)-V(X-ee n )>0, 

satisfies 

\Ah\ < CS 2 U n in B := Bjl(X*). 
Notice also that by Harnack inequality 

(4.25) |ifl<c f»xra. 

with C universal. Assume that 

g(X*)>V(X*). 
Then, in view of (j4~24l) and (|4~25)) 

h[X*) = g(X*) - V(X* - ee n ) > ceU n {X*). 
Hence by Harnack inequality, (|4.25j) and the condition CS 2 < e 

h > ceU n (X*) - C6 2 \\U n \\ L ~ {B) > c'eU n (X*) in Bj_(I*). 
Thus, using (|4.24|) we have that for r small enough 

h > c'esupK > V(X - (1 - r)ee n ) - V(X - ee n ) in 

from which our desired conclusion follows with any fj such that 2fj < min{?7/128, r}. 
Case 5. > ?y/4 and n P ^ 0. 

In this case we argue similarly as in the previous case but we need to make use 
of the boundary Harnack inequality. 

Assume that X* G {s > 0} and call X^ = (x*,0) the projection of X* onto 
{s = 0}. If e is small enough then it follows from (I4.24[) that the function 

h(X) := g(X) - V(X - ee„) > 0, 
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satisfies 

\Ah\ < CS 2 U n in B := Br (X *) n {s > 0}, 

for a universal constant C. Denote by Y* — Xq + j^e n and assume that 

g(Y*) > V(Y*). 

As in the previous case, by Harnack inequality 

(4.26) h>ceU n (Y*) in B^_(Y*). 

Now we argue similarly as in Lemma 14.41 
Denote by 

D := (B V/S (XS) \ B V/32 (Y*)) n {s > 0}. 
Let qi , q 2 satisfy in D 

A?i=0, Aq 2 = -1 
with boundary conditions respectively, 

q 1 = l ondB n/32 (Y*), q t = on d(B n/8 (X£) n {s > 0}) 

and 

q 2 = on 91?. 
By the maximum principle, in view of (I4.26[) we obtain that 
h > ceU n {Y*) qi - C5 2 q 2 in D. 

Moreover, 

q 1 > cq 2 in D D B v/16 (Xq). 

Hence using that CS 2 < e we get 

h(X) > c'eU n {Y*) qi (X) > ceU n (X) in B 7j/16 (X*) n {s > 0} 

where in the last inequality we used that U n {Y*)q\ is comparable to U n in view of 
boundary Harnack inequality. 

Now we use (|3.1|) and (|4.24[) to conclude 

= /i(x,x„ + i) > ce sup C/„(y,x n+ i) > ce sup 
B f (x *) s^(v *) 

> V(X - (1 - r)ee n ) - K(X - ee n ) in B^(X*) D B$_(X*). 

Then our desired statement holds for fj < min{r/2, 7y/64}. □ 

5. Improvement of flatness. 

In this section we prove our main Theorem 11.11 We start with the following 
quadratic improvement of flatness proposition. We show that if a solution g stays 
in a X 2+a neighborhood of a function 7 £ Vf in a ball B\ then in B v \, g is in a 
(Xf]) 2+a neighborhood of another function V in the same class. 

Proposition 5.1. Given a £ (0,1), there exist Xo,Vo & (0,1) and C > large 
depending on a and n, such that if g solves (|2.1[) . £ F(g) and g satisfies 

(5.1) V(X - \ 2+a e n ) < g{X) < V(X + X 2+a e n ), in B x with < A < A 

for V = VM.o,a,b £ Vi, then in a possibly different system of coordinates denoted by 
E = {ei, . . . , e„, e n+ i}, 

(5.2) V(X - (r/ A) 2+Q e„) < g(X) < V(X + (ry A) 2+Q e„), m B Vo \ 
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for some V — V§ a i (defined in (|3.3|) ) with S given in the E coordinates by 

S = {x n = i(x') T Mx'}, 

and 

||M-M||,|a-a|,|&-6| < CA a , a + b-trM = 0. 

Moreover, for any a £ (0, 1], the surfaces S andS separate in B a at most C(X a a 2 + 
X 1+a a). 

Proof. Let 770, C be the constants in Corollary 12.101 

The proof is by compactness. Assume that no such Ao exists, then we can find 
a sequence of A^'s, tending to 0, gf. and Vk satisfying (|5.1|) for which (|5.2[) fails. 
We rescale g^ and Vk ■ For simplicity of notation we drop the dependence on k and 
denote 

g x (X) = \- 1 / 2 g(\X), V X (X) = \- 1 / 2 V(\X), X e B x , 

Notice that 

Vx = V\M,0,\a,\b G V°, 



and 



Let 



Vx{X - X 1+a e n ) < g x (X) <Vx(X + \ 1+a e n ) in B x . 
e = A 1+Q , 5 = A 



and define 



(5.3) w x := 

e 

Thus by Proposition 



9\ ~ Vx Va -TVa / — ; , ^z^ 2 

»A = + = .9A e y. + 0( — 

e e ' A e 



and hence by Corollary 14 . 2 1 we get that wa converges uniformly to a Holder contin- 
uous function wq as k tends to 00 (and A — > 0), with wq(0) = and |u>o| < 1. 
We claim that wq is a viscosity solution of the linearized problem 



(5.4) 



{A(U n w ) = Q, in5 1/2 \P, 
jjV r i«o|=0, onB 1/2 nl. 



We start by showing that U n WQ is harmonic in B x / 2 \ P. 

Let (p be a smooth function which touches wq strictly from below at Xq € B X / 2 \P- 
We need to show that 

(5.5) A(U n <p)(X ) < 0. 

Since wx converges uniformly to wq in B1/2 \ P we conclude that there exist a 
sequence of constants cx — > and a sequence of points Xx <E -B1/2 \ P, X\ — > Xq 
such that ipx '■— £ (<^ + c a) + 7y touches gx by below at Xx for a sequence of A's 
tending to 0. 

Define the function ipx by the following identity 

(5.6) ^x(X-iix(X)e n )=U(X). 
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Then according to ()2.8[) ip\ touches g\ from below at Y\ = X\ — ip\(X\)e n E 
B^(g\). Thus, since g\ satisfies (|2 . 1 1) in B\ it follows that 

(5.7) A.tJj\(Y\) < 0. 

In a neighborhood of X , 7v A /A has bounded C k norms (depending on |Xo|) hence 
has also bounded C k norms. By Proposition [278] 

A^a = AA(C/„(^ a /A)) + 0(A 2 ) 

- A{U n jj x ) + 0(A 2 ) 

= A(U n (e^ + 1 ^))(X x ) + 0(X 2 ) 

= eA{U n $) + 0(A 2 ) 

where we have used that 

A(E/ n7v J = G. 

This can be checked either explicitly or by using Theorem 16.11 
In conclusion 

eA{U n ^){Y x ) + 0(X 2 ) <0. 
We divide by e = A 1+a and let A->0. Using that Y\ —> X we obtain 

A(U n tp)(X ) < 0, 

as desired. 

Next we need to show that 

\v r w \(x ) = o, x = 04 0,0) e B 1/2 nL, 

in the viscosity sense of Definition 12.91 

We argue by contradiction. Assume for simplicity (after a translation) that there 
exists a function <fi which touches wq by below at with 0(0) = and such that 

4>(X) = -x' + f3r + 0{\x'\ 2 + r 3 / 2 ), 

with 

(3>0. 

Then we can find constants er, f small and A large such that the polynomial 

q(X) =£-x'- ^\x'\ 2 + 2A{n - l)x n r 

touches 4> by below at in a tubular neighborhood N? = {\x'\ < f, r < f} of 0, 
with 

(j> - q > o > 0, on Nf \ N f/2 . 

This implies that 

(5.8) w - q> a > 0, onJV f \% 2 , 
and 

(5.9) w (0)-g(0) = 0. 

In particular, by continuity near the origin we can find a point X* such that 

(5.10) w (X*) -q(X*) < ^, X* e N f \ P close to 0. 

8 

Now, let us define 
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W X ■= VxM+Ael-e £' ,Xa,Xb+2eA(n~l) 

Then in view of Proposition 13.51 we have 

W x = eq + j Vx + 0(6 2 ) 
and moreover, W\ is a subsolution to our problem since e > i5 2 . 

Thus, from the uniform convergence of wx to w and (|5.8p we get that (for all A 
small) 

(5.11) 9x-W x =Wx ^ q + Q{ ^ ) > ° - m{Nf \ N )\ P , 

e e 2 

Similarly, from the uniform convergence of wx to wq and (|5.10l) we get that for k 
large 

(5.12) (gx-Wx)(X*) &tx * eNf \ P . 

e 4 

On the other hand, it follows from Lemma \2. 71 and (|5.1ip that 

gx-Wx a . 
> — m Nf \ ir. 

which contradicts (|5.12[) . 

In conclusion wo solves the linearized problem. Hence, by Corollary 12.101 since 
w (0) = 0, w satisfies 

(5.13) -~r^+ Q < Wo (X)-(t -x' + ±x /T M x'-^r 2 ~b x n r) < -^ 2+a in B Avo , 

for some rj 6 (0, 1) universal and with 

<z + 6 - irM = 0, |eo|,||Mo||,|ao|,|6o| <C. 
From the uniform convergence of wx to loo, we get that for all k large enough 

(5-14) -^< wx{ x)-^^<^r, 2 +a in B, no \P, 

with 

Tx '■= VxM—eMa,—e£o,Xa—eao,Xb—ebo ■ 

In conclusion, from the definition (15.31) of wx , we get 

(5-15) fx-^ 2 +a <~gx<fx + ^o +a , 

or 

Tx(X - e -vl +a e n ) < gx(X) < T X (X + | % 2+Q e„) in B 2vo . 
We rescale gx back from the ball B\ to Bx and obtain 

(5.16) T(X - y^ +Q e„) < g(X) < T(X + y^ +a e») in B 2Xm , 
with 

for 

S T := {z„ = \{x') T M T x' + £ T • a;'}. 
M T :=M-~M , £ T := -e£ , «t :=ffl-^o> &t := 6 - ^6 . 
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Next we show that in a different system of coordinates, called E, the function T 
can be approximated by VM T ,o,a T ,b T - 

Assume for simplicity that £t points in the e\ direction. Then we choose an 
orthogonal system of coordinates E := {ei, e-x, . . . , e n +i} with 

e, = e», ifi^l,n 

and e„ normal to St at 0. 

Notice that the E system of coordinates is obtained from the standard one after 
an orthogonal transformation of norm bounded by C|£y| which is smaller than Ce. 

A point in this system on coordinates is denoted by X. We let, 

S := {x n = i(x') T MTx'}, 

and we write S as a graph in the e n direction, that is 

S := {x n = h(x')}. 

We claim that in a ball of radius a the distance (in the e n direction) between St 
and S in B a is less that Cecr 2 , for any < a < 1. 

Indeed, since x = Ox with O orthogonal and \\0 — 1\\ < Ce, we obtain by implicit 
differentiation 

\\D 2 x ,h - M T \\ Lx{Bx) < Ce, V^/i(0) = fr. 
Thus in i?2r; A we have that the surfaces 

lie between 

S±eA % 2+Q e„ 

since Ce(?7oA) 2 < |A?7 2+Q! . 

In view of this inclusion, using that w ar ,fc T (t, s) is monotone in i, we obtain from 
(|5. 16|) the desired conclusion (|5.2[) with Af = Mt, a = a,T,b = br- 

Since the distance between St and £> in B a is less than C(— a 2 + ea) the proof 

A 

is finished. □ 

We can now prove our main Theorem ll.il In fact we show that under our flatness 
assumption, a solution g can be approximated in a C 2,a fashion by a function 
V g V° . 

Theorem 5.2. There exists e > small universal such that if g solves (|2.1j) m £?i 

(5.17) {x € £>i : x„ < -e} C {x G B 1 : g(x, 0) = 0} C {x G Bi : x n < e}, 

iften in an appropriate system of coordinates denoted by 

V(X - C\ 2+a e n ) < g{X) < V(X + CX 2+a e n ) in B x , for all < A < 1/C, 

for some V = VM ,o.a .b G w^/i C depending on n and a. In particular, 

F(g) n Z?i/2 is a C 2,a graph in the e n direction for any a G (0, 1). 

Proof. It suffices to prove the theorem for any fixed a G (0, 1) for some e(a), C(a) 
depending on a. The dependence of e on a can be easily removed by fixing e := e(a), 
say with ci = 1/2. Then by the conclusion (|5.2j) for a, appropriate rescalings of g 
satisfy the flatness assumption ()5.17|) also for e(a) for any a G (0, 1). 
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By Lemma EU the rescaling 

g li (X) = ii-igQiX) 

satisfies 

U(X - r 2+a e n ) < 9il {X) < U(X + T 2+a e n ) in B x , 

provided that e, fj, are chosen small depending on t < Ao, with Ao the universal 
constant in Proposition [lO] and r small universal to be made precise later. Thus 
satisfies in B T the hypotheses of Proposition 15.11 with M = 0, a = 0, b = 0. Then 
we can apply Proposition 15 . II repeatedly for all := ttjq since by choosing r small 
enough we can guarantee that J^k ^ T k — 1 anc ^ hence the corresponding Mfe, ak, 
bk have always norm less than 1. Thus we obtain 

(5.18) V Sk , ak , bk (X - r 2+a e k n ) < g^X) < V Sk , ak , bk (X + r 2+a e k n ) in B Tk . 

Using that Sk and Sk+i separate (in the e„-direction) in B„ at most C(r k a 2 + 
T k v+a u) we conclude that as k — > oo, the paraboloids Sk converge uniformly in B\ 
to a limit parabolid <S*. Moreover, 5* also separates from Sk in B a by at most 
C{r k a 2 + T^ +a a) in the e* direction where e* is the normal to S* at the origin. 
Finally, as k — > oo, ak a*, bk — > £>*, with 

- a*|, - 6, | < Ct£. 

Now notice that in B2 Tk , the paraboloids Sk and S» separate at most Cr 2+a , 
thus we can apply Lemma 13.71 and use the inequality (|5.18l) to obtain 

Vs„ at ,bAX - Cr 2 k +a e* n ) < 9lM (X) < V s ,^j } ,(X + Crl +a e* n ), in B Tk . 

Rescaling back we obtain the desired claim. 



□ 



6. The regularity of the linearized problem 
We recall that the linearized problem associated to (I2.ip is 



(6.1) 
where 



A(U n h) = 0, in B x \ P, 
\V r h\ = 0, on B x n L, 



IV7 u/v "\ t h(xQ,x n ,s) h(x o ,0,0) 2 2 , 
v r /i A := hm - , r = x„ + s 

(x„,s)->(0,0) r 



2 



In this section we obtain a second order expansion near the origin for a solution 
h to ED- 



Theorem 6.1. Let h be a solution to (|6.1[) swc/i f/ia£ < 1. Then h satisfies 

(6.2) \h(X) - (h(0) +do-x' + \{x') T M Q x' - ^r 2 - 6 rz n )| < C\X\\ 

for some ao,bo,^o, Mq with |£q|, |oo|j I&o|j ll^oll < C, C universal and 

<iq + bo — trMo = 0. 
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Proof. This proof is a refinement of Theorem 8.1 in [DRj where the authors obtained 
a first order expansion for h, in particular 

(6.3) \h{X)-h(X )\ =O(\X-X \), X eL. 

Also in |DR| it is shown that h and its derivatives of all orders in the x' direction 
are Holder continuous with norm controlled by a universal constant in -B1/2 (see 
Corollary 8.7.) 

We wish to prove that 

(6.4) \h(x',x n ,x) - h(x' ,0,0) + -y^f 2 + b(x')x n r\ < CY 3 , (x',0,0) € B 1/2 DL, 

with C universal and h{-, 0, 0), a, b smooth functions of x' . 
The function h solves 

A(U n h) = in Bi \ P, 
and since U n is independent on x' we can rewrite this equation as 

(6.5) ^x n , s (U n h) = -U n A x ,h. 

Moreover, since A x >h solves the same linear problem as h then any estimate for h 
also holds for A x >h. 

For each fixed x', we investigate the 2-dimensional problem 

A(U t h) = U t f, in B 1/2 \{i<0,s = 0}cl 2 

with h, f € C°'P. Without loss of generality, for a fixed x' we may assume 
h(x', 0, 0) = 0. Thus in view of (|6.3p . the function 

H := U t h 

is continuous in B1/2 C K 2 and satisfies 

AH = U t .f in B-l/2 \ {t < 0,s = 0}, H = on B 1/2 n {t < 0, s = 0}. 
Now, we consider the holomorphic transformation z — > ^z 2 

$:(C,y)^(i, s ) = (i(C 2 - y 2 ),Cy) 

which maps B\ n {C > 0} into B 1 / 2 \ {t < 0, s = 0} and call 

ft(C,!/) = M*.«). f(C,y) = f(t,s), H(C,y) = H(t,a) 
with (f, 0), the polar coordinates in the (£, y) plane. Then, easy computations show 
that 

(6.6) AH=Cf inS 1 n{C>0}, H((,y) = ^h, 
and 

H = on {( = 0}. 

Since the right-hand side is in C '^ and ft,, / have the same regularity, we conclude 
from repeatedly applying Lemma ROI below that h, f G C°° with 

ll/llo*./9(B+ a )' V4c^{B+ 2 ) ^ c ( k >P)- 

Notice that we can reflect H oddly and h, f evenly across {( — 0} and the resulting 
functions will still solve (I6.6[) in B\. Moreover from our assumptions, / and h are 
even with respect to y. Thus, we conclude that the Taylor polynomials for /, h 
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around the origin, are polynomials in Q 2 , y 2 . Now we use the Taylor expansion for 
H around 0, which is odd with respect to C and even in y, that is 

H((, y) = ((d Q + d^ 2 + d 2 y 2 + 0(f 4 )) 

with 

6di + 2da = /(0) = /(0). 

Thus, 

HC, v) = C 2 (do + diC 2 + d 2 y 2 + 0(f 4 )). 
In terms of the (£, s) coordinates this means that 

6 6 
h[t, s) = 2r{d + 2d!r(cos -) 2 + 2d 2 r(sin -) 2 ) + 0(r 3 ) = 



= 2r(d + {d x + d 2 )r + {d x - d 2 )t) + 0(r 3 ) 

In conclusion, 



2d r - |r 2 - btr + 0(r 3 ). 



\h{X) - h(x', 0, 0) - 2d (x')r + ^p. r 2 + b(x')x n r\ < Cr 3 in B 1/2 
with C universal, 

IMU°°({|x'|<l/2}) ; II^IIl~({|x'|<1/2}) ^ C > 

and 

a + b = A x ,h(x',0,0). 
Since h solves (|6.1j) we must also have do(x') — and hence 

\h(X) - h(x', 0, 0) + + b(x')x n r\ < Cr 3 , in B 1/2 . 

Notice that a, b are smooth functions of x' with all order derivatives bounded 
by appropriate universal constants. Indeed due to the linearity of the problem it 
is easy to see that D^,a,D^,,b are the corresponding a and b for D^,h. Writing 
the Taylor expansions at for h{x' , 0, 0) up to order 2 and a, b up to order 1 with 
a = a(0), b = 6(0) we get 



h(X) - (h(0) +£o-x'+ ^(x'fMox 1 - yr 2 - b x n r^j 



< C\X\ 3 . 

□ 

In our proof above we used the following easy lemma. 



Lemma 6.2. Let H = H(C, y) be a function defined on B^ C R 2 , which vanishes 
continuously on {( = 0}. If H e C k - a (B+), k £ N, a G (0,1], then C^H E 
Ck-h»(B+), and 

nr^iic— < \\h\\c^. 

Proof. Since H(0, y) — we see that 

C 1 H((,y)= f H c (t(,y)dt 
Jo 

and the lemma follows easily by taking derivatives in the equality above. □ 
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7. Basic properties of a solution g. 

We collect here some useful general facts about solutions g to our free boundary 
problem (|2.1[) . such as C 1 ' 2 -optimal regularity, asymptotic expansion near regular 
points of the free boundary and compactness. 

First we recall some notation. Let v £ C(Bi) be a non- negative function. We 
denote by 

B+(v) := B 1 \ {(x, 0) : v(x, 0) = 0} C R" +1 

and by 

F(v) :=flt»(B^(w)nBi)nBi Cl" 
Also, we denote by P the half-hyperplane 

P :={X £ R n+1 :x n <0,s = 0}. 

Given a C 2 surface S in R™ -1 , we often work with functions of the form V = 
Vs,a,6 (see Definition [33]). We remark that we can still apply the boundary Harnack 
inequality with V in a neighborhood of S since in this set V is comparable with a 
harmonic function H with F(H) = S. 

Indeed, after a dilation we may assume that V = Vs.a.b £ V#, that is the curva- 
tures of S in B 2 and \a\, \b\ are bounded by S small, universal. Let 

Vl := Vs,a-2n8A, V% '■= Vs, a +2nS,b 

and notice that V\ is a supersolution and Vi is a subsolution in B\ (see Proposition 
[32| • Also 

1/2V2 < V < 2Vi < 2V 2 , 

hence there exists H between 1/2 V2 and 2Vi, with 1/4F < H < 41/, harmonic 
in {ff > 0} and F(#) = 5. 

We obtain the following version of the boundary Harnack inequality. 

Lemma 7.1. Let V := Vs, a ,b G Vs a > f or some small Sq universal and with £ 
S. Let w £ C{B\) be a non-negative function which is harmonic in B^(w). If 
B+(V) C B+(w) then 

w > cw(~e n )V, in B 1/2 . 

IfB+(V) C Bf (w) then 

w <C\\w\\ L ^ (Bl) V, in B 1/2 . 

Proof. Let w be the harmonic function in B^ 4 (V) with boundary value w on dB^/^ 
and w = where {V = 0}. 

If B^(V) C B^(w) then in view of the observation above we can apply the 
boundary Harnack inequality with V and conclude that 

w > w > cw(^e n )V, in B 1/2 . 

On the other hand, 

w{— e n ) > c inf w. 

2 S 1/4 (fe„)naB 3 /4 

Using that w and w coincide on dB 3 / 4 together with Harnack inequality we obtain 
our desired estimate. 
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If B?(V) C BJ{w) then 

w < w < Cw(-e n )V, in B 1/2 . 

On the other hand, 

w(-e n ) < \\w\\ L oa {B3/i) = \\w\\ L oo {B3/i) 
which yields our conclusion. 

□ 

An immediate consequence is the following useful lemma. 

Lemma 7.2. Let V = Vs a .b G Va be a subsolution in B\, for some small universal 
So and with £ S. If w is harmonic in B~^(w) and B^(V) C B^(w) and 

w > V — e in B\ , 

then 

w > (1 - Ce)^ m S 1/2 . 

Proof. Let g : -> 1 be the harmonic function in B3/4 n -B^F) which has 

boundary values q = 1 on dB 3 / 4 and q = on the set where = 0. From our 
hypotheses on w and the maximum principle we obtain 

w > V — eq in -B3/4. 

On the other hand by Lemma [7.11 since B^, 4 (V) = B^ 4 (q) we have q < CV in 
B1/2, which together with the inequality above implies the desired result. 

□ 

Remark 7.3. From the proof of Lemma [7721 we see that if the hypotheses on w hold 
only outside of the ball Si/g, i.e 

w > V — e on Bi\B 1 / s , w harmonic in B+(V)\B 1 / S 

then the conclusion holds in the shell -B3/4 \ -Bi/4- 

Next we prove optimal C 1 / 2 regularity for viscosity solutions. 

Lemma 7.4 (C 1 ^ 2 -Optimal regularity). Assume g solves (|2.ip in B\ andO £ F(g)- 
Then 

g(x,0)<C\d(x)\ 1 / 2 mB 1/2 
where d(x) represents the distance from x to F(g). Also 

Us llcV2( Bl/2 ) < C(l + g(~e n+x )). 

Proof. The first assertion follows in a standard way from the free boundary con- 
dition. By scaling, we need to show that if g is defined in B2, € F(g) and 
Bi{e n ) C B+(g) then u(e n ) < C for some large C universal. 

By a rescaled version of Lemma 17.11 and Harnack inequality we have that in a 
neighborhood of 0, 

.9 > cg(e n )V s ,2n,o, S = dBi(e„) 
with Vs,2n,o a subsolution near 0. The free boundary condition gives 1 > cg{e n ) 
which provides a bound for g(e n ). 
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For the second inequality we write 

9 = 9o +5i in D : = B 3/4 H {s > 0}, 
with go,gi harmonic in D and satisfying the following boundary conditions 
go = g on {s = 0} n dD, g a = on {s > 0} H &D, 

c/i = on {s = 0}n<9L>, ,9i =.g on {s > 0} n dD. 
From our estimate for g on {s — 0}, we obtain 

ll.9o||ci/2(B 1/2 nD) < C||3llci/2(g 3/4 ) < C, 
which together with the bound 

1 1 

IMIcv=>(B 1/3 nr>) < Cffi(2 e «+i) ^ CffC 2^+!)' 

gives the desired conclusion. 



□ 



Next we prove that if F(g) admits a tangent ball at either from the positive or 
from the zero phase, then g has an asymptotic expansion of order od^l 1 / 2 ). This 
expansion also justifies our definition of viscosity solution to the free boundary 
problem (|2.ip . We remark however that this expansion holds also for an arbitrary 
harmonic function w which does not necessarily satisfy the free boundary condition. 

Lemma 7.5 (Expansion at regular points from one side). Let w € C X I 2 (B\) be 
1/2-Holder continuous, w > 0, with w harmonic in B^(w). If 

€ F(w), B 1/2 (l/2e n ) cB+(w), 

then 

w = all + od^l 1 / 2 ), for some a > 0. 
The same conclusion holds for some a > if 

B 1/2 (-l/2e n ) C {w = 0}. 

Proof. We define 

a := inf liminf — (tv). 
First we notice that a > 0. Indeed, by a rescaled version of Lemma l7.ll 

w > cw(-e n )V s , q,o, S = dBi(^e n ) 

near the origin, for some c > 0. This implies that a > cw{ — e n ) > 0. 

Assume by contradiction that the conclusion of the lemma does not hold with 
this choice of a. Then there exist 5i > and a sequence of points y& — ^ such that 

(7.1) \w(y k )-aU(y k )\ > Jil^ 1 ' 2 . 

Since w is 1/2-Holder continuous on B\, the rescalings 

w k {x) := \yk\~*w(\Vk\x), 
are uniformly 1/2 Holder continuous and after passing to a subsequence we can 
assume that w k converge uniformly on compact sets to a limiting function w* € 
C(R n ). We obtain 

w* > all, Aw» =0 in K™ \ P, 
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and in view of ()7.1[) there exists a point y* , | y* | = 1 such that 

u>*(j/*) > aU(y*) + 5\. 
Using boundary Harnack inequality we find 

(7.2) w* > a(l + <5 2 )t/ in Si, 
for some ^ > small. Now we let 

^■7,0,52,0 

and we notice that V is subharmonic in B\ (by Proposition [372]) and satisfies 
V(X)=v S2 ,o(t,s) = (l + jp)U(t,s) 

<(l+ 5 j-)U(t,s)<(l+^)U(x n ,s). 

Thus ||Z2]) gives, 

(7.3) w»>a(l + ^)U inBi. 

From the existence of a tangent ball at the origin included in {w > 0} we see that 
for all large k, Wk is harmonic in the set where {V > 0}. Thus we conclude from 
([73]) that in Bi 

Wk > a(l + ~tW ~ £fe > f° r somc e k — > 0. 
By Lemma I7T21 we find that for all large k, 

Wk > (1 - C— )a(l + -i)V > a(l + -§-)V in B 1/2 . 

a 4 8 ' 

This implies that for any v ^ P 

liminf —(tv) = liminf — -(tv) > a(l H — -), 

t^o+ U t->o+ U 8 

which contradicts the minimality of a. □ 

Remark 7.6. If we assume that F(w) admits a uniform tangent ball from its side 
at all points in B1/2 then the hypothesis w € C x ^ 2 (B\j^) is satisfied and therefore 
w has an expansion at all points in F(w) nBi/4. Indeed, by Lemma l7.ll we know 
that 

w < C\\w\\ L o,V 9Br(xo)fifi 
with dB r (xo) a tangent sphere to F(w) from the side, and this implies 

w(x) < C\\w\\ L oo dist{x, {w = 0}) 1/2 , \/x e Bi/4, 

which gives w £ C 1 ^ 2 (B 1 / 4 ). 

In general, the term od^l 1 / 2 ) in the expansion for w can be improved in o(U) 
in the non-tangential direction to F(w). For example assume that £ F(w) £ C 2 
and e n is the normal to F(w) at which points towards the positive phase. Then 
the non-tangential limit 

r w 
lim — = a 

x£C,x~>0 U 

where C C M" \ P is a cone whose closure does not contain L = {x n = 0, s = 0}. 
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Indeed, by Lemma \7. 5 1 and Remark 17.61 we have that w = a U + odX] 1 / 2 ). Now 
the limit above follows by applying boundary Harnack inequality for U and w in 
the sets C\ n (B r \ B r / 2 ) for all r small, where 

d :={\x'\ >fj,\{x n ,a)\} 

is such that C C C\ U {0}. 

Remark 7.7. In the definition of viscosity solutions for our free boundary problem 
(see Definition 2.3) we can restrict the test functions only to the class of subsolutions 
and supersolutions of the form cVs, a ,b- 

Precisely we say that g is a solution to (|2.1[) if 

1) A 5 = 0in B+(g); 

2) for any point X £ F(g) there exists no Vs, a ,b such that in a neighborhood of 
Xq, Vs : a,b is a subsolution and 

g > aVs, a ,bi for some a > 1 

with S touching strictly F(g) at Xq from the positive side. 
Analogously there is no supersolution Vs, a ,b such that 

g < aVs ; a,b-> f° r some a < 1 

and S touches strictly F(g) at X from the side. 

In order to prove this statement we need to show that if we can touch g by below 
at a point Xq £ F(g) with a comparison subsolution v as in Definition 2.2, then we 
can touch also with a subsolution aVs t a,b as above. A similar statement holds for 
supersolutions. 

Assume for simplicity that Xq = 0, e n is normal to F(v) at and g > v in B\. 
Let v be the harmonic replacement for v in B±(v). In view of of Remark 1 7. 6 1 

v = all + o(\X\ 1/2 ), for some a > 1. 

We claim that a > 1. Indeed, v — v > is superharmonic in B^(v) and vanishes 
continuously on {v = 0}ni?i. If v — v = 0, then our claim follows from the definition 
of a comparison subsolution. Otherwise, by the boundary Harnack inequality 

v — v > av 

in a neighborhood of the origin, for some a > 0. Thus v > v/(l — a) near the origin 
and again the claim follows from the expansion of v at the origin. 

— 1/2 

The rescalings Vf- = r k vyrf-x) converge uniformly on compact sets to all, with 
a > 1 . As in the proof of Lemma 17.51 we obtain that there exists 6 small such that 
for all large k, 

v k >v~ (i + <y)v&/,o,*,o 

and F(V) touches strictly F(vk) at the origin from the positive side. Rescaling 
back we obtain the desired conclusion. 

Next we prove a compactness result for viscosity solutions to (|2.ip whose free 
boundaries converge in the Hausdorff distance. 

Proposition 7.8 (Compactness). Assume gu solve (|2.ip and converge uniformly 
to g* in B\, and {gu = 0} converges in the Hausdorff distance to {<?* = 0}. Then 
<?* solves (|2.1[) as well. 
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Proof. Clearly is harmonic in (<?*). In view of Remark l7.7l we need to check say 
that if £ F(g*) there exists no subsolution VM.o,a,b such that in a neighborhood 
of 0, 

.g* > aVM,o, a .b, for some a > 1 

and F(V) touches strictly F(g*) at from the positive side. A similar statement 
can be checked also for supersolutions. 

Assume by contradiction that such &V — VM,o,a,b exists. Then after a dilation 
we may assume that V € Vs for some small S and V is a subsolution in B\ . 

For any e > there exists a > such that for all \t\ < a and all large fc's 

W t (X) :=aV(X + te n ) <g k ~e in B 1 , 

and 

F(W- a ) C B+(g k ), F(W t ) \ B 1/s C B+{g k ). 
By Lemma 17.21 and Remark 17.31 we obtain that 

9k > (1 - Ce)W- a in B 1/2 

and 

9k > (1 - Ce)W t in B 3/4 \ B 1/4 , for all |t| < a. 

By choosing e small (depending on a) we see that the functions Wt are strict 
subsolutions to our free boundary problem, and hence the inequality above can be 
extended in the interior (see Lemma T2.7|) i.e., 

9k>{l-Ce)W t inB 1/2 . 

Writing this for t = a we see that {g k = 0} stays outside a neighborhood of the 
origin and we contradict the convergence in the Hausdorff distance to {<?* =0}. 

□ 

We conclude this section by showing that our flatness assumption on the free 
boundary F(g), implies closeness of g and U. 

Lemma 7.9. Assume g solves (|2.1|) . Given any S > there exist e > and fi > 
depending on S such that if 

(7.4) {i£8i:x„< -e} C {x G B x : g{x, 0) = 0}c{ie6i: x n < e}, 
then 

U(X-fiSe n )<g(X)<U(X + fj,6e n ) mB^. 

Proof. The proof is by compactness. Assume by contradiction that a sequence of 
functions g k satisfies the hypotheses with e k —> but the conclusion does not hold. 

Notice that by Harnack inequality g k {e n +i/2) is bounded be a multiple of g k {e n /2) 
which in view of Lemma 1 7. 41 is bounded by a universal constant. Hence by the sec- 
ond claim in Lemma T7. 41 the g k s have uniformly bounded C 1 / 2 norms on compact 
subsets of B\. After passing to a subsequence we can assume that g k converges 
uniformly on compact sets of B\ to a function with 

(7.5) A 5 * = infli\P, 5* = onPnBi. 

By Remark I7.6[ is C 1 / 2 . Moreover, the derivatives of g* in the x' direction 
satisfy again (|7.5[) and we obtain 

\\D^g*\\c^ { B 1/2) <c(fi). 
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Now we can separate the variables and write 

^x n ,s9* — 

and we can argue as in Theorem 16. II to obtain 

\g*{X)~aU{X)\<C\X\^ 

with C universal. 

We now want to apply Proposition l7.8l to conclude that solves (|2.1[) and hence 
a = 1. To do so, we must guarantee that > in B\ \ P. Otherwise = and 
hence ||gfc||L°°(s 1/2 ) — > 0. Let B k ■— Bi/$(xk) be a ball tangent to F{g k ) from the 
zero side at some point y k € -Bi/8- Then, since ||3fc||L oo (Si /2 ) ~~ ^ 0; we have by 
Lemma 17.11 

9k < VkVdB k ,o,0> witn °k -> 0. 
This contradicts the free boundary condition for g k at y k . 
In conclusion, g* solves (|2.ip and 

\ 9 4X)-u(x)\<c\x\V 2 

with C universal. 
Rescaling we find 

\g k ^(X)-U(X)\<Cv in£? 2 , with g k ,„{X):=iT l l 2 g k [jjX). 

Thus 

9kA X ) > U ( X ) ~Cn> U(X - e feM - 1 e„) - C/i. 
Now we use that F(g k u) C {|^ n | < an d obtain by Lemma [7T21 that in B x 

9k,„ > (1 - C/i)C/(X - efc/i-^) > C/(X - (e^- 1 + C M K) 

where the last inequality follows once more from boundary Harnack inequality. A 
similar inequality bounds g k ^ by above. We choose \i small depending on 8 and 
obtain that g k ^ satisfies the conclusion of the theorem 

U(X - 5e n ) < ^ 2 g k {^X) < U{X + 5e n ), 
and we reach a contradiction. □ 
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